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In this paper we consider a spatially fiat Friedmann-Robertson-Wallcer (FRW) cosmological model 
whit cosmological constant, containing a stiff ffuid and a classical Dirac field. The proposed cos- 
mological scenario describes the evolution of effective dark matter and dark energy components 
reproducing, with the help of that effective multifluid configuration, the quintessential behavior. 
We find the value of the scale factor where the effective dark energy component crosses the phantom 
divide. The model we introduce, which can be considered as a modified ACDM one, is characterized 
by a set of parameters which may be constrained by the astrophysical observations available up to 
date. 
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I. INTRODUCTION 

According to the standard cosmology the total en- 
ergy density of the Universe is dominated today by both 
dark matter and dark energy densities. The dark mat- 
ter, which includes all components with nearly vanishing 
pressure, has an attractive gravitational effect like usual 
pressureless matter and neither absorbs nor emit radia- 
tion. The dark energy component in general is consid- 
ered as a kind of vacuum energy with negative pressure 
and is homogeneously distributed and, unlike dark mat- 
ter, is not concentrated in the galactic halos and in the 
clusters of galaxies. The observational data provide com- 
pelling evidence for the existence of dark energy which 
dominates the present day Universe and accelerates its 
expansion. 

In principle, any matter content which violates the 
strong energy condition and possesses a positive energy 
density and negative pressure, may cause the dark en- 
ergy effect of repulsive gravitation. So the main problem 
of the modern cosmology is to identify this form of dark 
energy that dominates the universe today. 

In the literature the most popular candidates are cos- 
mological constant A, quintessence and phantom mat- 
ter. Their equation of state is given hy w = p/ p, where 
w = —1, w > ~l and w < —I respectively. Dark en- 
ergy composed of just a cosmological term A is fully 
consistent with existing observational data. However, 
these data do not exclude the possibility of explaining 
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the observed acceleration with the help of phantom mat- 
ter. The cosmological constant can be associated with a 
time independent dark energy density, the energy density 
of quintessence scales down with the cosmic expansion, 
and the energy density of phantom matter increases with 
the expansion of the universe. 

Mostly, the attention has been paid to dark energy as 
high energy scalar fields, characterized by a time varying 
equation of state, for which the potential of the scalar 
field plays an important role. Among scalar field models 
we can enumerate quintessence models ll| , Chameleon 
fields [3| , K-essence Q jChaplygin gases [j| , tachyons [1] , 
phantom dark energy |6|, etc. 

In general the crossing of the phantom divide cannot 
be achieved with a unique scalar field 7] this fact has 
motivated a lot of activity oriented toward different ways 
to realize such crossing [8|. For instance, in Ref. Bit 
was explored the so called kinetic k-essence models 
i.e. cosmological models with several k-fields in which 
the Lagrangian does not depend on the fields themselves 
but only on their derivatives. It was shown that the dark 
energy equation of state transits from a conventional to 
a phantom type matter. Note that formally, one can 
get the phantom matter with the help of a scalar field by 
switching the sign of kinetic energy of the standard scalar 
field Lagrangian @. So that the energy density pph = 
-(l/2)$2 + ]/($) and the pressure = -(l/2)$2 _ 
of the phantom field leads to p^h +Pph ~ < 0, 
violating the weak energy condition. 

In the Universe nearly 70% of the energy is in the form 
of dark energy. Baryonic matter amounts to only 3 — 4%, 
while the rest of the matter (27 %) is believed to be in 
the form of a non-luminous component of non-baryonic 
nature with a dust like equation of state {w = 0) known 
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as cold dark matter (CDM). In this case, if the dark 
energy is composed just by a cosmological constant, then 
this scenario is called ACDM model. 

Below, we analyze a FRW universe having cosmolog- 
ical constant and filled with a stiff fluid and a classical 
Dirac field (CDF). With this matter configuration, we 
will see that the FRW universe evolves from a non ac- 
celerated stage at early times to an accelerated scenario 
at late times recovering the standard ACDM cosmology. 
The CDF may be justified by an important property: in 
a spatially flat homogeneous and isotropic FRW space- 
time it behaves as a " perfect fluid" with a energy density, 
not necessarily, positive definite. This pressuless "per- 
fect fluid" can be seen as a kind of " dust" . In particular 
motivated by the fact that the dark matter is generally 
modelled as a system of collisionless particles [ill, [i3 1 
have the possibility of giving to cold dark matter content 
an origin based on the nature of the CDF. On the other 
hand, the stiff fluid is an important component because at 
early times, it could describe the shear dominated phase 
of a possible initial anisotropic scenario, dominating upon 
the remaining components of the model. 

The organization of the paper is as follows: In Sec. 
II we present the dynamical field equations for a FRW 
cosmological model with a matter source composed by a 
stiff fluid and a CDF. In Sec. Ill the behavior of the dark 
energy component is studied. In Sec. IV we conclude 
with some remarks. 



of the flat space-time constant Dirac matrices 7* as 

v/3 
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where the Dirac matrices 7* can be written with the Pauli 
matrices as 
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(6) 



,(7) 



with / the identity 2x2 matrix. The symbol Vi — 
di + Yii denotes the spinorial covariant derivatives, being 
the spinorial connection defined by ViF^ = 0. Then 
it leads to 



So = 0, 



/3- 



(8) 



The constant a will be associated, later on, with the total 
observable matter. 

Coming back to the Dirac equation it takes the 
form 



r° ( a, + 



m = 0. 



(9) 



II. DYNAMICAL FIELD EQUATIONS 

We shall adopt a spatially flat, homogeneous and 
isotropic spacetime described by the FRW metric 



(1) 



where ait) is the scale factor. The spacetime contains a 
cosmic fluid composed by (i) a stiff fluid ps — Pso / and 
(ii) a homogeneous classical Dirac field tp. The Einstein- 
Dirac equations are 



3i72 _ A = ^ 



PD, 



(2) 
(3) 



where H = d/a is the Hubble expansion rate, a is a 
constant and the dot denote differentiation with respect 
to the cosmological time. Here pu represents the energy 
density of the CDF. 

The dynamical equation for the CDF in curved space- 
time can be obtained using the vicrbein formalism. So, 
F* are the generalized Dirac matrices, which satisfy the 
anticommutation relations 



(4) 



with the metrics tensor 5*'^' and / the identity 4x4 matrix. 
They can be defined in terms of the usual representation 



Restricting ourselves to the metric ([T|), the general so- 
lution of the latter equation consistent with Eq. ^ is 
given by 



62 

die'"'* 



(10) 



with arbitrary complex coefficients bi, 62, c?i and d2. The 
only nonvanishing component of the energy-momentum 
tensor for the CDF is 



TS, = ^{\h\' + \b2\'-\d,\'-\d,\')^P^, (11) 
where poo Oi{b'^ - d^), 5^ = + |52p and = 



Ml 



For positive values of poo this source for- 



mally behaves as a perfect fluid representing a classical 
dust. However, the CDF will allow us to extend the anal- 
ysis for negative values of the energy density. Here, we 
restrict us to the physical sector d^ < 5^ and without loss 
of generality we choose b^ — 1; this means that po > 0- 
From Eq. pT|) we see that the energy density of the 
CDF is given by po — poo/a.^- Thus the Eq. ([2]) can be 
rewritten in the following form 



Px 



(12) 
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FIG. 1: We show the behavior of the scale factor (straight 
Una) and its derivatives a (dotted line) and a (dashed line). 
From this it is clear that this cosmological scenario exhibits 
an accelerated expansion since there is a stage where a > 0. 



FIG. 2: We show the behavior of as a function of time. 
The dashed line represents the limit case = am, i.e. a^d* = 
4Apso , while straight and dotted lines represent a typical case 
satisfying the condition Oc > Om since a^d* < 4Apso. 



where 



a 



Px 



PsO 



A. 



(13) 
(14) 



We shall associate the positive part of the Dirac energy- 
momentum tensor with the pressureless matter giving 
rise to the total ("true") observable matter p^, while 
one may assume that its negative part, along with the 
stiff fluid and the cosmological constant constitute the 
effective dark energy component p^- 

The general solution of the Einstein equation (fT2)) with 
sources and ^1)) takes the form 



a(l -d^) 




(15) 



where we have setting the initial singularity at t — 0. In 
the Fig [Hit is showed the behavior of the scale factor and 
its derivatives. 



III. DARK ENERGY EVOLUTION 

In order to have a positive px, we choose the parame- 
ters of the models according to the following restriction 



(16) 



Since a{t) is an increasing function of time, this effective 
dark energy component decreases up to it reaches a min- 
imum value pxc = A — a^d^/ipso at Uc = {2pso/ acPY^'^ 
where the dark component crosses the phantom divide 
and begins to increase with time (see Fig. [2]). Funda- 
mentally, the effective dark energy component crosses the 
phantom divide due to presence of the c?-parameter. In 
fact, the cosmological evolution of dark energy depends 
on the negative term — ad^/a"^, see Eq. (fT4)) . because it 
produces a minimum at Pxio-c) — 0, showing the impor- 
tance of considering the CDF as a source of the Einstein 
equation. 

Assuming that total matter and dark energy are cou- 
pled only gravitationally, then they are conserved sepa- 
rately, so we have that 



p„i + 3Hpm = 0, 
Px + 3i7(l + Wx)px = 0, 



(17) 
(18) 



where we have assumed the equation of state Px — ^xPx 
for the dark energy component. Taking into account that 
the dark energy component has a variable state parame- 
ter Wx = Wx{a), we define as — (pso/A)^/^ the value 
of the scale factor where this equation of state coincides 
with the matter one, that is Wx = 0. So that, the above 
restriction a^d* < 4Apso now becomes Um < cic- In terms 



FIG. 3: We show the behavior of the dark energy state pa- 
rameter Wx (dashed hne) and the state parameter for the full 
source content (straight line). It is clear that the dark 
energy violates the dominant energy condition while the full 
source content does not violate it. Note also that the dark 
energy component remains in the phantom region as it enters 
to it after crossing the phantom divide. 



of the above parameters, Oc and a^, the dark energy state 
parameter Wx can be written as 



1 - {a/a„if 



1 - 2{a/acf + (a/a„)6 



(19) 



On the other hand we can consider the state equation for 
the fuU source content. For this 2-component system we 
define the total pressure = w^p^, where = Px and 
the total energy density Pj. — p,n + Px- This implies that 

(20) 



1 + r' 

where r = Pm/ Px- In the Fig. Owe compare the behaviors 
of the state parameters of the dark energy and the full 
source content. 

Let us consider more in detail the ratio of energy den- 
sities of matter and dark energy. The defined above ratio 
takes the form 



(21) 



Aa^ — ad^a^ + psa 



It is easy to show that this ratio has a maximum at a = 
flc = (Pso/A)^^*^ (at this point the (fr/da^ < 0). This 
maximum value is given by 



2pso - ad'^ v^Pso/A 



(22) 



FIG. 4: We show the behavior of the ratio r = pm/px as a 
function of the scale factor a. The dashed line represents the 
ACDM cosmology, where r oc The solid line represents 

the case where r > 1 and so at the beginning dark energy 
dominates, then at the range a £ (a-,a+) the dark matter 
component dominates over dark energy, and for a > a+ the 
dark energy component dominates again. The doted line rep- 
resents a scenario where the dark energy always dominates 
over the dark matter, since rmax < 1. 



It is clear that for cosmological scenarios where r^nax < 1 
the dark energy component always dominates over the 
dark matter during all cosmological evolution. Thus in 
order to have stages where dark matter dominates over 
dark energy we have to require that Vmax > 1- So in this 
case we would have two values for the scale factor where 
the energy density of dark matter equals the energy den- 
sity of the dark energy: 



a± 



a{l + d^) ± v/a2(i + ^2)2_4A^^g 



2A 



1/3 



(23) 



Thus for cosmological scenarios where r^ax > 1 at the 
beginning the dark energy dominates over the dark mat- 
ter until the scale factor reaches the value a = a_ where 
the dark matter energy density equals the energy density 
of dark energy and it begins to dominate. This stage of 
domination of the dark matter is prolonged to the mo- 
ment where the scale factor reaches the value a — and 
the dark energy starts to dominate again over the dark 
matter (see Fig. S]). 
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A. Constraints on cosmological parameters 

The proposed scenario is characterized by four param- 
eters which may be constrained by the astrophysical ob- 
servations available up to date. Since we have considered 
flat FRW cosmological scenarios the dimensionless den- 
sity parameters are constrained today as 



^m,0 + ^x.O — 1- 

From Eqs. dnj-llll) we have that 

2 _ a a(P 



(24) 



3H' = — + ^ 
a3 flS 



A. 



3 , ... (25) 

Evaluating it today (where we set a = 1) we have that 

Pent = 3Hq = a + Psq - a(f + A, (26) 
so the two dimensionless density parameters are given by 
Praia = 1) _ a 



^m,Q — 



Pcrit 



+ PsO ~ cd? + A ' 



(27) 



n^,= Pli^^^El^^l^_ (28) 



Pcrit 



a + PsQ — adP + A 



Now it may be shown that in general, for a flat FRW 
cosmology deceleration parameter q is given by 



9 = 



aa 1 3 
^ " 2 ^ 2^ 



(29) 



Taking into account that we have = px = uJxPx and 
Pt — Pm + Px the deceleration parameter is given by 



1 3 

9=2 + - ^m), 



(30) 



and evaluating it today (i.e. a — \) and using p9p we 
obtain 

«(1 - <fi){l - 2(7o) + 2/9,0(2 - go) = 2A(1 + go)- (31) 

Thus, for accelerated scenarios, q < 0, we require a posi- 
tive cosmological constant. Other constraint may be in- 
troduced by taking into account the moment when the 
universe has started to accelerate again. In other words 
this is related to the moment when the Universe starts 
violating the strong energy condition, i.e. p + p>Q and 
/o + 3p > 0. So we must require the inequality p + ip < 0. 
Now from condition a = and the equivalent Friedmann 
equation 



a 1 , 
-^-7:{p + 3p), 
a 6 



(32) 



we conclude that p + 3p = 0, which implies that p^^ + 
Px + 'iuJxPx = 0, obtaining the condition 

a(l - + Zaccf + 4p,o(l + Zaccf = 2A, (33) 



where the equation 1/a = (1 -I- z) was used. Here Zacc 
is the value of the redshift where the universe starts to 
accelerate again. 

In conclusion we have the four conditions (P7|) . (|28p . 
([3T|) and p3|) for the four parameters a, p^Oj d and A of 
our model. 

Note that from Eqs. (|3ip and ([55)1 we have that 



Pdo = a(l - d^) = Kpso, 



(34) 



where 



(1 - 2qo) - (l + qo)z 

Since poo is related to the energy density of the CDF we 
must to require that K > 0. So from Eqs. (HH), (pS)) . 
(ISTI) and (l33l) we have that 



PsO = 



2K 



A 



an.o[^(2 + z) + 2(l + 2z2)]' 



m^zjK + 4z) 
K{2 + z) + 2{l + 2z^)'' 



(36) 
(37) 

(38) 
(39) 



where we have used the Eqs. ((24)) and ((34)) . 

Now the four model parameters need to be constrained. 
We made this by using the Eqs. ((M)) - ((55)) and by con- 
sidering the increasing bulk of observational data that 
have been accumulated in the past decade. The present 
expansion rate of the universe is measured by the Hub- 
ble constant. From the final results of the Hubble Space 
Telescope Key Project [l^ to measure the Hubble con- 
stant we know that its present value is constrained to 
be i7o = 72 ± 8 kms^^ Mpc"\ or equivalently H^^ = 
9.776/i^^ Gyr, where /i is a dimensionless quantity and 
0.64 <h<0.8 12]. 

Now assuming a flat universe, i.e. Eq. ((24)) is valid, 
Perlmutter et al. found that the dimensionless den- 
sity parameter fim.o may be constrained to be 0.3, 
implying that Qx,o 0.7 [H, [lH|, and the present 
day deceleration parameter qo may be constrained to be 
-1< go < -0.64 

For consistency we need also to compare the age of 
the Universe determined from our model with the age of 
the oldest stellar populations, requiring that the Universe 
is older than these stellar populations. Specifically, the 
age of the universe to is constrained to be fo > 11-12 
Gyr 

So let us calculate the age of the universe from Fried- 
mann equation ((25p . This equation may be written as 



PsO + aa^{l - (f) + a* 
a + PsO ~ ad? + A 



(40) 
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obtaining Eq. ([24]) when the expression (|40)) is evaluated 
today. Then the age of the universe may be written as 

^0 - Jo H{l+z) - 

HE Jo V «(i-rf^)(i+2)^+p,o(i+2)«+A iT+^y- 

In the Table U we include some values obtained from 
Eqs. ((36)) ~ ((39)) for the parameters a, d, pso, A and X cor- 
responding to some given values of the parameters Hq, 
qo and Zacc- For the Hubble parameter Hq is consid- 
ered both possible values i?o- for h — 0.64 and i?o+ for 
h = 0.8, and for the matter dimensionless density param- 
eter we have taken rim.o = 0.3. The last two columns rep- 
resent the age of the universe determined from the model 
parameters for h = 0.8 and h = 0.64 respectively. Clearly 
in the proposed model there are configurations which are 



allowed from oldest stellar age since there exist combi- 
nations of the parameters a, d, Pso, A which give > 
11-12 Gyr satisfying the stellar population constraints. 



Ho 


qo 


^acc 


Q (in units of SHq ) 




Pso (in units of SHq) 


A (in units of 3-ffo) 


K 


Ho+,to (Gyr) 


Ho-, to (Gyr) 


Ho+ 


-0.68 


0.587 


3.86 X 10"^ 


0.38 


1.7 X 10"" 


10"^ 


14 


10 


12 


Ho- 


-0.68 


0.587 


2.85 X 10"^ 


0.38 


1.3 X 10"'' 


7.6 X 10"^ 


14 


10 


12 


Ho+ 


-0.68 


0.94 


3.86 X 10"^ 


0.29 


1.05 X 10"^ 


10~^ 


2619 


12 


15.7 


Ho- 


-0.68 


0.94 


2.85 X 10"^ 


0.29 


7.73 X 10"^ 


7.5 X 10"^ 


2619 


12.5 


15.7 


Ho+ 


-0.9 


0.6 


3.9 X 10"^ 


0.964 


3.5 X 10"'' 


1.2 X 10"^ 


0.38 


10 


12.7 


Ho- 


-0.9 


0.6 


2.8 X 10"^ 


0.964 


2.6 X 10"" 


9.1 X 10"^ 


0.38 


10 


12.7 


Ho+ 


-0.9 


1 


3.9 X 10"^ 


0.815 


7.2 X 10"^ 


1.2 X 10~^ 


9.9 


12.5 


15.6 


Ho- 


-0.9 


1 


2.8 X 10"^ 


0.815 


5.3 X 10"^ 


8.9 X 10"^ 


9.9 


12.5 


15.6 



TABLE I: In this table we show some values of the model 
parameters obtained for given Ho, qo and Zacc {^m,o = 0.3). 



B. The effect of the d— parameter 

Now it is interesting to get some insights concerning 
the nature of the d-parameter for studying the effect on 
the cosmological evolution of the negative term —ad^/a^ 
in the Eq. for the dark energy component px ■ It can 
be shown that in the proposed cosmological scenario the 
dominant energy condition (DEC) is violated thanks to 
the presence of this parameter. Effectively, if d = the 
dark energy state parameter Wx and the state parameter 
of the full source content are given by 

"^x = r-^J W„ — — T, (42) 

respectively. From these expressions we see that always 
— 1 < Wa; < 1 which implies that now the dark energy 
component satisfies DEC, as well as w^. Their general 
behavior is showed in the Fig. O (compare with Fig. [3]). 

However, the fulfilment of the DEC does not imply that 
the universe has a decelerated expansion. From Eq. (|32p 



we may write that 

— = -{Pt + ^Vt) = 

and putting d = we see that the accelerated expansion 
is realized if a > aacc = ((« + + 32Apso)/4A)i/3. 
Other property of the d-parameter to be considered is its 
effect on the deceleration parameter go- From Eq. (|38p . 
which is independent of the Hubble parameter Hq, and 
using Eqs. (j35p we can express the deceleration parameter 
qq as a function of the Zacc obtaining 

1 — z 

qoiZacc) = -2 _^ {6d'^n,n,0 Z 

-6 /2,„^o z + 4z + 3d2/2™,o - 3 /?™,o -I- 4) , (44) 

where as before z = (1 -I- Zacc)^- We can see that in this 
case qo{zacc) rapidly tends to the value 

go,oo = -l + |(l-d2)0„,o, (45) 
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FIG. 5: We show the behavior of Wx and ui^, as a function of 
the scale factor. The dashed hne represents the behavior of 
the dark energy state parameter, while straight line represents 
the behavior of the state parameter of the full source content. 
Both satisfy the DEC since in this case d = Q. 



FIG. 6: We show the behavior of the deceleration parameter 
qo as a function of Zacc for three values of the parameter d 
(0, 1/2, 1) with Slm.o = 0.3. We can see that the deceleration 
parameter rapidly tends to the values (—0.55, —0.6625, —1) 
respectively. 



obtaining for d — and rim,o — 0.3 the value go^oo = 
—0.55, and from this value the deceleration parameter 
reaches the value (70,00 = — 1 for d « l(see Fig. [6]). So 
the parameter d affects directly the range of validity of 
the deceleration parameter which is constrained to be in 
the range — 1 < (70 < 0. Note that the value (70,00 = — 1 
for (i = 1 is independent of the value of the dimensionless 
density parameter 0.m,Q- 



IV. CONCLUDING REMARKS 

Since today the observations constraint the value of 
w to be close to w = — 1, we have considered broader 
cosmological scenarios in which the equation of state of 
dark energy changes with time. The two principal in- 
gredients of the model are a stiff fluid which dominates 
at early time and a CDF. The positive part of the lat- 
ter was associated with a dark matter component while 
its negative part was considered as a part of the dark 
energy component and was the responsible that the ef- 
fective dark energy density crossed the phantom divide 
(see and compare the Figs. |3]and[5]). At the end, this 
cosmological model becomes accelerated recovering the 
standard ACDM cosmology. 

In general, the model may be seen as a continuation of 
the inflation era. In the inflationary paradigm the scalar 
field driven by the potential V{^), generates the in- 



flationary stage. In the slow roll limit 6^ << V^($), with 
Lo^ K, —1, we have a super luminal expansion while in 
the kinetic-energy dominated limit $^ >> V{^), with 
LJ<^ « 1, we have a stiff matter scenario characterized by 
a subluminal expansion. Taking into account that our 
model has a variable equation of state, we can think it as 
a transient model which interpolates smoothly between 
different barotropic eras, as for instance, radiation dom- 
inated era, matter dominated era an so on. In other 
words, from Eqs. (fT9l) - ((2T|) we see that — > 1 (and 
^ 1) for a ^ implying that the energy density px 
behaves like 1 / and matching, after inflation, with the 
kinetic-energy mode of the scalar field /?$ cx 1/a®. Now 
from Eq. ^9]) we see that passes trough a radiation 
dominated stage (i.e. = 1/3) for 



^rad 



22/3 




(46) 



behaving like px ~ V^-rad ^^'^ dominating over After 
that, at a = Om, we have Wx — and the effective dark 
component behaves as a pressureless source obtaining a 
matter-dominated stage. Finally the model evolves from 
this state to a vacuum-energy dominated scenario. It 
is interesting to note that in the a matter-dominated 
stage if the condition (fTB|) is fulfilled then the total energy 
density is given by = a{l — d^)/a^-^ -\- 2A, and if am = 
ttc then p^ = a/a^, being px — 0, implying that we have 
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at this stage only the dark matter component. 

The above results indicate that a cosmological sce- 
nario based on a CDF component and the effective multi- 
fluid configuration can, in certain cases, reproduce the 
quintessential behavior (see Figs. [3] and [5])). In fact, the 
state parameter of the total matter content — 1 < < I 
is constrained the same as the state parameter of the 
scalar field in quintessence models. In this manner we 
avoid the use of scalar fields and particular classes of po- 
tentials for describing the dark energy component. 

Finally, all the parameters of the model has been ex- 
pressed in terms of the observable quantities which may 
be constrained by the astrophysical observational data. 
In effect, in the Table U some values of the model param- 
eters a, d, pso and A were included, which correspond 
to some given values of the parameters Hq, qo, Zacc and 
f^m,o constrained by astrophysical observations. 
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